Coincidence of the oscillations in the dipole transition and in the persistent current of 

narrow quantum rings with two electrons 
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The fractional Aharonov-Bohm oscillation (FABO) of narrow quantum rings with two electrons 
has been studied and has been explained in an analytical way, the evolution of the period and 
amplitudes against the magnetic field can be exactly described. Furthermore, the dipole transition 
of the ground state was found to have essentially two frequencies, their difference appears as an 
oscillation matching the oscillation of the persistent current exactly. A number of equalities relating 
the observables and dynamical parameters have been found. 
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Quantum rings containing only a few electrons can 
be now fabricated in laboratories^'^. When a magnetic 
field B is applied, interesting physical phenomena, e.g., 
Aharonov-Bohm oscillation (ABO) and fractional ABO 
(FABO)of the ground state (GS) energy Eo and persis- 
tent current Jo, have been observed 2-4, -j-j^g -j-j^g_ 
oretical aspect, a number of calculations based on exact 
diagonalization^"®, local-spin-density approximation^'^^, 
and the diffusion Monte Carlo method^^ have been 
performed. These calculations can in general repro- 
duce the experimental data. For examples, in the 
calculation of 4-clcctron ring^'^^, the period of oscilla- 
tion $0/4 found in experiments was recovered (<i>o — 
hc/eisthefluxquantum) . 

In addition to the oscillations in Eo and Jo, the oscil- 
lation in the optical properties is noticeable. ^^^7^ jjj ^j-^jg 
paper a new kind of oscillation found in the dipole tran- 
sition of two-electron (2-e) narrow rings is reported. The 
emitted (absorbed) photon of the dipole transition of the 
GS was found to have essentially two energies, their dif- 
ference is exactly equal to hJo, where h is the Planck's 
constant. In other words the difference of the two photon 
energies appears as an oscillation which matches exactly 
the oscillation of Jo- This finding is approved by both 
numerical calculation and analytical analysis as follows. 

The narrow 2-e ring is first considered as one- 
dimensional, then the effect of the width of the ring is 
further evaluated afterward. The Hamiltonian reads 
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where m* the effective mass, 9j the azimuthal an- 
gle of the j — th electron, $ — nR'^B/^o, where B 
is a magnetic field perpendicular to the plane of the 
ring, V12 the e-e Coulomb interaction, Hzeeman = 
—Sz^J'^ the well known Zeeman energy where Sz is 
the Z-component of the total spin S, and fi ~ ^^i^ 1 



where g* is the effective g- factor and /is is the Bohr 
magneton. The interaction is adjusted as ^ Vyi = 

e^/{2e^(F + R'^ sin^{{ei - 92)12))-^ , where e is the di- 
electric constant and the parameter d is introduced to 
account for the effect of finite thickness of the ring. 

We first perform a numerical calculation so that all 
related quantities can be evaluated quantitatively, m* = 
0.063me, e = 12.4 (for InGaAs), d = 0.05i? , and the 
units meV , nm , Tesla and <I>o are used. Accordingly, 
G = 604.8/i?^and ^l^33.53/R^. 

A set of basis functions (j)k,k.2 = e'e^i^i+'^^^^VSTr is in- 
troduced to diagonalize the Hamiltonian, where fci and 
/c2 must be integers to assure the periodicity, the sum of 
ki and k2 is just the total orbital angular momentum L. 
4>kxk2 must be further (anti-)symmetrizcd when S — 0(1). 
When about three thousand basis functions are adopted, 
accurate solutions (at least six effective digits) can be ob- 
tained. The low- lying spectrum is plotted in Fig.l, where 
the oscillation of the GS energy and the transition of the 
GS angular momentum Lo can be clearly seen. 

Let Oc = {O2 + ei)/2 , and 6I2 - 6*1. Then 
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where Hcoii = ^G{-i-g^ + ^^)--\-iize 
2G{—i-^Y + ^12, they are for the collective and internal 
motions, respectively. Our numerical results lead to the 
following points. 

(i) Separability: The separability of one-dimensional 
ring is well known^. However, for the convenience of the 
following description, it is briefiy summarized as follows. 
Each eigenenergy E can be exactly divided as a sum of 
three terms 
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E = \G{L -f 2$)2 + Em - 
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where the first term is the kinetic energy of collective 
motion, Ei^t is the internal energy . 

Since the basis functions can be rewritten as 
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= e*^^c7e*5(fe-fei)v'/27r (4) 

the spatial part of each eigenstate ^ is strictly separa- 
ble as = -^e*^ ^^"^Hnt where the first part describes 
the collective motion, while ^int is a normalized internal 
state depending only on ip. In particular, both Eint and 
V'int do not depend on B (or $). 




FIG. 1: Low-lying levels of a 2-e ring against $/3>o in the 
FABO region. When $/ is positive, Lo is negative, the num- 
bers by the curves are —Lo- 
in) Classification oftpint- When L is even (odd), (^2 — 
ki)/2 is an integer (half- integer), thus the period of ip as 
shown in (4) is 27r (47r). Therefore, the periodicity of the 
internal states have two choices. In fact, the difference in 
the periodicity is closely related to the dependence of the 
domains of the new variables 9c and (p, this point has 
been discussed in detail in ref.[14,15]. Let Q = (— 1)^, 
then the four cases (Q, S) = (1,0), (-1,0), (-1,1) and (1,1) 
are associated with four types of states labeled by a, b, c, 
and d , respectively. The internal states of Type a 
arc denoted as tpa, ipa' , • • • a-ud the associated internal 
energies as Ea < Ea*, ■ ■ ■ and so on. Examples of 
V'mt and Eint are plotted in Fig. 2 and listed in Table 1, 
respectively. 

Table 1, The lowest and second lowest internal 
energies (in meV) of Type a to d, R = 30nm. 



Type 


a 


b 


c 


d 


Eint 


2.626 


4.247 


2.630 


4.272 


^int 


6.342 


8.912 


6.435 


9.158 



Due to the c-c repulsion, a dumbbell shape (DB), i.e., 
= 180°, is advantageous in energy because the two 
electrons are farther away from each other meanwhile. 
However, a rotation of this geometry by tt is equivalent 
to an interchange of particles, these operations will cre- 
ate the factors (—1)^ and (— l)*^, respectively, from the 
wave function. Therefore, the equivalence leads to a con- 
straint, accordingly the DB is allowed only for the states 
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FIG. 2: Four types of ipint against ip , R = 4Qnm. The lowest 
three of each type are shown, the higher state has more nodes. 

with L + S even, i.e., only for Type a and c. Other- 
wise, the states would have an inherent node at the DB 
and therefore be higher in energy as shown in Table 1, 
where Ea « Eb, E^ « Ed, and Ea » Ec- In Fig.2 
the patterns of Type a arc one-to-one similar to Type c 
, they all have a peak at the DB. On the contrary, all 
those of Type (6) and(rf) have the inherent node at the 
DB. It is noticeable that Type b and c are not continuous 
at = and 27r due to their periods are not equal to 27r. 
It was found that the internal states of all the GSs are 
either tpa or tpc without exceptions because the favorable 
DB is allowed in them. When the dynamical parame- 
ters vary in reasonable ranges, the qualitative features of 
Fig.2 remain the same. 

According to (3), an appropriate Lq would be chosen to 
minimize the GS energy. When $ increases, Lo will un- 
dergo even-odd transitions repeatedly and become more 
negative as shown in Fig.l. Correspondingly, the total 
spin So undergoes singlet-triplet transitions, and ipa and 
ijjc appear in the GS alternatively. However, due to the 
Zeeman effect, when $ is larger than a critical value ^crit 
, only 5o = 1 states will be dominant, and accordingly 
only V'c will appear in the GS. The region $ < (>) ^crit 
is caUed the FABO (ABO) region. 

(iii) Persistent Current: Let J\ be the current of the 
particle ei. The expression of Ji is well known."'' However, 
since it does not depend on the azimuthal angle, it equals 
to its average over 6i . Thus the total current J = J1 + J2 
is 

J=J^gJd0^d92 [**(-zJ^-zJj + 2<l>)* + c.c.] (5) 

where g = h/{m*R^). Using the arguments 9c and (p 
and making use of the separability, the integration over 
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6c and <^ can be performed. Thus we have 

J = g{L + 2$)/27r (6) 

This equation demonstrates exphcitly the mechanism 
of the oscillation of the persistent current, it is caused 
by the step-by-stcp transition of L during the increase 
of $. Examples of J are shown in Fig. 3, where each 
stronger oscillation (associated with a L odd and S = 1 
GS) is followed by a weaker oscillation (associated with 
a L even and S = Q GS). 
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FIG. 3: The oscillation of the persistent current and the two 
photon energies of the ground states against $/$o in the 
FABO region. Tiio unit of current is lO^ 'C/R, where C 
is the velocity of light. In the lowest panel, the black square 
(white circle) denotes hujj^ {hui-), namely, the energy associ- 
ated with Lo to I/O + 1 {Lo — 1) transition. 

(iv) Relations among the internal states: Define 

Cm = e""(''i-'''^)+ e^^iO-2-<}c) = 2cos(mv?/2). By an- 
alyzing the numerical data, we found 

A^(OiV'a) = V'6 + and N{Oii)c) = V'd + (7) 

where N is the operator of normalization, both and 
are very small functions and depend on the dynamical 
parameters very weakly. E.g., when R varies from 30 to 
90, the weights of and vary from 0.0004 to 0.0002. 
They arc so small that in fact can be; neglected. Since 
Ox contains a node at the DB, it must cause a change of 
type from a to 6, or from c to d. Thus it is not surprising 
that (7) holds. Since Oi is the operator of the dipolc 
transition (see below), eq.(7) provides an additional rule 
of selection as discussed later. 

(v) Dipole transition: The probability of dipole tran- 
sition reads P^^) ± = ^ {uj±/c)^R'^ l^(oj±l^' where u!± 
is the frequency of the photon, 



=SL,,„L,^,±i{4il^m4:^,) (8) 

where (/) and (o) denote the final and initial states, 
respectively, the signs ± are associated with L(j) = 
Lie) ±1- 

Let the initial state be the GS with Lg, then i/'l^l 
must be ipa or depending on Lq is even or odd . 
Let a denotes the type of the initial state. Due to (7) 

, {4il^\Oi\4:l) = hf),<^ < Oi^alOiVa where 
S(f)^a implies that the final state must be ipi, (ipd) if a = a 
(c), otherwise the amplitude is zero. Thus, due to the ad- 
ditional rule of selection eq.(7), the dipole strength of the 
GS is completely concentrated in two final states having 
iy(j) = iy(o) ± 1 and both having the same internal state 
specified by eq.(7). Accordingly, only the photons with 
the two energies 

fko± = E(^f)±. - E(^g) 

= G[i(l±2(L„ + 2$)) + A,/G] (9) 

can be emitted (absorbed), where = — Ea or 
Ed — Ec depending on a = a or (c). The oscillation 
of huj± is plotted in the lowest panel of Fig. 3. It turns 
out that Aa/G depends on R very weakly, thus Sw± is 
nearly proportional to Accordingly, a smaller ring 

will have a larger probability of transition with a higher 
energy. 

(vi) FABO region: The oscillation in this region is com- 
plicated as shown in Fig.l and 3. It is noted that the GS 
energy (3), persistent current (6), and the photon ener- 
gies (9) all contain the factor Lo + 2$, thus their FABO 
arc completely in phase and have the same mechanism 
caused by the transition of Lo against In Fig.f the 
abscissa $ can be divided into segments, in each the 
GS has a specific Lo and the GS energy is given by a 
piece of a parabolic curve. The segment is called an 
even (odd) segment if Lo is even (odd). At the border of 
two neighboring segments the two GS energies arc equal. 
From the equality and based on (3), the right and left 
boundaries of the segment with Lo can be obtained as 

<^right{Lo) = (l-(/z/2G)2)-l[ 1-I^{E,-Ea)/G^- 

2Lo+(-l)^<'(2(^e-£;a)+M(io-l/2))/G]/4 (10) 

<i>leftiLo) = (1 - (Ai/2G)2)-lhl - fi{E, - Ea)/G^ - 

2Lo - (-l)^°(2(i?, - K) + Ai(io + l/2))/G]/4 (11) 

where Lo < and ^right{Lo) = ^ieft{Lo - 1), arises 
from the Hzeeman- The length of the segment reads 

dL„ = ^r^ght{Lo) " ^left{Lo) = (1 - (m/2G)2)-1[1 + 

{-l)^^{2{E,-Ea)+^Lo)IG]l2 (12) 
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which is related to the period of the FABO. When $ 
increases, the magnitude of Lo would increase. Since (iLo 
is negative, it is clear from eq.(12) that the length of even 
(odd) segments would become shorter (longer) when $ 
increases. 

The location of a segment with a given Lg can be 

known from the inequality ^leftiLo) < $ < ^rightiLo)- 
Once the relation between and the segments of $ is 
clear, every details of the FABO can be analytically and 
exactly explained via the cq.(3), (6), and (9). In particu- 
lar, the extrema in each segment can be known by giving 
^ = bright OT ^left- For an example, the maximal cur- 
rent is g{Lo + 2^right)/'2TT. Incidentally, the minimum of 
the GS energy in a segment is i^min = E^— ji^ /SG+iJ,Lo/2 
(if 5o = 1), or just equal to Ea (if 5o = 0). 

It is noted that Ec - Ea (cf. Table 1) and /i/G (it 
is 0.0554 in our case) are both small. When <l> is small 
the magnitude of \Lo\ would be also small . In this case 
eq.(12) leads to cIl^ « 1/2, i.e., the period is a half of 
the one of the normal ABO. In fact, (12) provides an 
quantitative description of the variation of the period of 
the FABO. 

(vii) ABO region: When $ becomes sufficiently large, 
Lo will become very negative, the even segments will 
disappear due to their lengths dj^^ < . We can de- 
fine a critical odd integer Lcrit so that cIl^^h-i < 
while ^Lerit-M > 0' thereby the critical flux separating 
the FABO and ABO region can be deflned as 



^crit = ^left{Lcrit) 



(13) 



Once $ > ^crit, Lo remains odd and the system keeps 
polarized. Let Ix be the largest even integer smaller 
than -{G + 2{Ec - Ea))/n. It turns out from eq.(12) 
that Lcrit = -fx + 1- With our parameters, Lcru = — 19 
and accordingly ^cHt = 9.003 ( refer to Fig.l). Both 
Lcrit and '^crit depend on R very weakly, but sensitively 
on the effective mass m* . 

In the ABO region ($ > $crit), eqs.(lO) to (12) do not 
hold. Instead we have bright = —{Lo — l)/2, ^left = 
-(Lo -I- l)/2, and = 1. Thus the normal ABO re- 
covers. Evaluated from (6), the magnitude of current is 
from — g/27r to g/2TT (for a comparison, it is from —g/Air 
to g/AiT for 1-e rings). From (9) the photon energies 
is from Ac — G/2 to Ac + 3G/2, at the same time fuv- 
is from Ac + 3G/2 to Ac - G/2. 

(viii) Relations between the photon energies and other 
physical quantities: Due to (7), the emitted (absorbed) 
dipole photon has only two frequencies , therefore it is 
meaningful to define A^;^ = fi.(w+ — lv-). Directly from 
(9) and (6), we have 



hJn 



(14) 



where h is the Planck's constant and Jo is the persis- 
tent current of the GS. To compare with 1-e rings, the 



latter has An^ = 2/iJoEq.(14) demonstrates that the os- 
cillation of Afi^ and the oscillation of Jo are matched 
with each other exactly, they keep strictly proportional 
to each other during the variation of <I>. 

The maxima of Aroj measured in the ABO and FABO 
regions, respectively, read 



(Aft,) 



AB 



2G 



(15) 



(16) 



Obviously, (15) provides a way to determine G, m*can 
be thereby obtained. (16) can be rewritten as 

Ec-E, = (G-mLc)/2-(2G-m)/(4G)(A,J^4^ (17) 
This equation can be used to determine Ec — Ea- Fur- 
thermore, we define 



Tnu> = h{uj+ +uj-) = G + 2A„ 



(18) 



Once G has been known, (18) can be used to determine 
Eb—Ea and Ed—Ec- Since the spectrum can be generated 
from the internal energies via (3), the evolutions of the 
spectrum and the persistent current against $ can be 
understood simply by measuring the photon energies. 
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FIG. 4: Evolution of /i Jo (solid line) and A^uj (dotted line) 
against B for a 2-e ring with ra = 50 and rt = 120nm 

(ix) Effect of the width: We now consider a two- 
dimensional model in which the two electrons are strictly 
confined in an annular region by a potential U{r), which 
is zero \i Va < r < Vb or is infinite otherwise. Under 
this model wc have performed numerical calculation to 
obtain A^ui and hJo, where Jo is now the total angular 
current inside the ring (from Va to rfe). The result is 
shown in Fig. 4 where ra = 50 and = 120 arc assumed, 
and the two quantities are slightly different from each 
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other. However, when the width becomes smaller, say 
Tb — fa < 30, the two curves overlap. Thus (14) works 
not only for one-dimensional but also for two-dimensional 
narrow rings. Let us define f = h/ -J m* ( A/ij^)^^^. 
For one-dimensional rings and from (15), we have f= R, 
where R is the radius of the ring. For two-dimensional 
rings, it was found from our numerical calculation that 
r {rb + ra)/2 if n - ra < 30. E.g., when n = 100 
and ra = 70, r =85.03. When n = 100 and Va = 90, r 
=95.00. Thus (15) works also well for two-dimensional 
narrow rings if the i? in G is replaced by the average 
radius. 

It is noted that the band-structure and related optical 
properties of 2-e rings have already been studied in de- 
tail by Wendler and coauthors^^. They classify the eigen- 
states according to their radial motion, relative angular 
motion, and collective rotation. In our paper the relative 
angular motion is further classified into four types ac- 
cording to the inherent nodal structures and periodicity 
of their wave functions, i.e., according to whether the DB 
shape is allowed and whether the wave function is contin- 
uous at if = 27r. The DB-accessibility turns out to be im- 
portant because it affects the eigenenergies decisively. In 
fact, the classification of states based on inherent nodal 
structures was found to be crucial in atomic physics, 
this would be also true in two-dimensional systems. Fur- 
thermore, the rule of selection for the dipole transition 
has been proposed in ref. [18]. In our paper, an addi- 
tional rule (namely,eq.(7)) is further proposed based on 
the possible transition of internal structures. This rule 
would affect the dipole spectrum seriously because the 
emission (absorption) is thereby concentrated into two 
frequencies. The difference of these two frequencies turns 
out to be proportional to the persistent current. There- 
fore the measurement of this difference can be used to 
determine the magnitude of the current. 

In summary, we have studied the FABO both analyt- 
ically and numerically. The analytical formalism pro- 
vides not only a base for qualitative understanding, but 
also provides a number of formulae for quantitative de- 
scription. The domain of <I> is divided into segments, 
each corresponds to a Lq. This division describes ex- 
actly how Lo would transit against \Phi, which causes 
directly the FABO. Thereby the variation of the period 
and amplitude of the oscillation of the GS energy, persis- 
tent current, and the frequencies of dipole transition in 
the FABO region can be described exactly. A number of 
equalities to relate the physical quantities and dynamical 
parameters have been found. In particular, a new oscilla- 
tion, namely, the oscillation of A^;^ was found to match 
exactly the oscillation of Jo- Since the photon energies 
can be more accurately measured, other observables and 



parameters can be thereby determined via the equali- 
ties. Since the separability of the Hamiltonian and the 
existence of inherent nodes are common, the above de- 
scription can be more or less generalized to iV— electron 
rings, this deserves to be further studied. 
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